Interplay of spin density wave and superconductivity with different pairing symmetry 
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A model study for the coexistence of the spin density wave and superconductivity is presented. 
With reference to the recent angle resolved photo emmission experimental data in high Tc cuprates, 
presence of the nested pieces of bands is assumed. The single band Hubbard model, therefore, 
when treated within the Hatree-Fock mean field theory leads to a spin density wave (SDW) ground 
state. The superconductivity (SC) is assumed to be due to a generalised attractive potential with a 
separable form without specifying to any particular origin. It therefore allows a comparative study 
of the coexistence of superconductivity of different order parameter symmetry with the spin density 
wave state. We find that the phase diagram, comprising of the amplitudes of the respective gaps 
(SC and SDW) Vs. band filling resembles to that of the high Tc cuprates only when the order 
parameter of the superconducting phase has d-wave symmetry. Thermal variation of different order 
parameters (e.g, SC and SDW) also show interesting coexistence and reentrance behaviors that are 
consistent with experimental observations, specially for the borocarbides. 
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I. INTRODUCTION 

The coexistence of magnetism and superconductivity 
has intrigued both experimentahsts and theoreticians for 
a long time. It is well known that the superconductivity 
(SC) and spin density wave (SDW) co-exist in a large 
classes of systems e.g., organic layered superconductors 
§, the Bechgaard salts {TMTSF)X, with X = PFq, 
AsFq, CIO4 etc. 1^, heavy fermion systems Q, high 
temperature superconductors Q and recently discovered 
borocarbides 

In case of the organic superconductors belonging to the 
family {TMTSF)2X {TMTSF = tetramethyltetrase- 
lena fulvalence and X = PFq, AsFq etc.) show coexis- 
tence of superconductivity and spin density wave like an- 
tiferromagnetism (SDW) at low temperature in moderate 
pressures (~ 7 Kbar) |6). Many compounds {X — PFq, 
AsFq, SBFq etc.) exhibit metal - insulator transition at 
12 - 17 K under ambient pressure which is caused by 
spin density wave (SDW) transition. This transition is 
suppressed by moderate pressures that results in super- 
conductivity. However, the nature of this coexistence is 
very much different from what would encounter in case of 
ternary compounds. In organic superconductors, there is 
only one conduction band and the antifcrromagnetic or- 
dering destroys the Fermi surface (FS) due to its 2kF 
periodicity which results in an insulating phase. In high- 
Tc cuprates also due to the low dimensionality aspects 
there exists similar periodicity due to nesting of the FS 
0, and leads to FS instability by forming the SDW state. 
On the other hand, in intermetallic compounds the mag- 



netism is due to localized / electrons and superconduc- 
tivity arises due to the conduction electrons. The or- 
ganic superconductors in turn, are essentially quasi one 
dimensional compounds with weak interaction between 
the adjacent chains. Our main concern in the present 
paper will however be to discuss the coexistence of an- 
tiferromagnetism (SDW) and superconductivity (SC) in 
the cuprates. 

In the high temperature cuprate superconductors 
La2-x(Sr, Ba)xCu04_y and YBa2Cu307_y strong anti- 
ferromagnetic spin order with a commensurate wave vec- 
tor Q = (7r/a, 7r/b) has been observed in neutron Q and 
Raman scattering |^] experiments. Moreover, in contrast 
to the conventional low Tc superconductors, the coher- 
ence length of high- Tc materials is short (10 to 20 pA) ; 
the average magnetic field seen by the cooper pairs need 
not be zero. Thus it appears that the antiferromagnetism 
and superconductivity are intimately related in these ma- 
terials, in contrast to the magnetism opposing supercon- 
ductivity in conventional superconductors. 

In the context of the influence of spin correlation on su- 
perconductivity in oxide systems there exists two limiting 
points of view ; (i) localized picture where U ^ W and 
(m) itinerant picture U <W ; where U is the strength of 
the onsite coulomb repulsion and W is the band width. 
In the large U limit, the valence electrons are localized on 
each Cu-site. To lowest order in t/U there is an antifcr- 
romagnetic exchange interaction J of the order of /U , 
which produces spin ordering [t being the hopping inte- 
gral) [|lO| . In the second point of view, U is assumed to be 
small so that the Mott Hubbard transition does not takes 
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place and Bloch states form a convenient basis to describe 
the system. Now for a perfect square lattice and at half- 
filling, the bands get nested with certain nesting wave 
vector Q that determines a new super periodicity lead- 
ing to spin density wave ground state which may or may 
not be comensurate with the underlying lattice. The spin 
density wave ground state causes an energy gap 21!^$ dw 
at the Fermi surface of the undoped material, leading to 
an insulator. Here we shall confine our selves to the latter 
point of view (i.e. U <^ W). 

The SDW ground state is obtained from the single 
band Hubbard model within the Hatree-Fock approxi- 
mation and assuming that the nesting of the FS exists 
only in certain directions of the FS. The directions of the 
FS where nesting exists will be instable with respect to 
the SDW formation whereas the superconducting insta- 
bility may occur in the rest of the part of the FS, pro- 
vided there exists some attractive interaction between 
the quasi-particles mediated by some boson exchange. 
Now, as regards the high- Tc cuprates are concerned the 
knowledge of the exact mechanism for the high tempera- 
ture superconductivity is yet to be known. In confirmity 
with the role of electronic correlations, the parent (un- 
doped) cuprate materials are always antiferromagnetic 
(SDW state in our case) but the long range AF order 
vanishes sharply on doping the system. However, there is 
some evidence that the AF fluctuations exists even in the 
superconducting phase. Hence there have been propos- 
als among the others that the pairing might as well arise 
due to the exchange of quanta of AF fluctuations [ pT| . 
In that case, the superconducting order parameter will 
have a d-wave symmetry. On the other hand, the SC in 
the system may even be caused by the quanta of the fluc- 
tuations of the collective (amplitude and phase) modes 
of the SDW state, present in our model, as proposed by 
Schrieffer et al., and others orelse the SC may 
even be due to normal phonon mediated mechanism. In 
either of these cases the nature of the superconducting 
order parameter symmetry could be different depending 
on the nature of the origin of the pairing mechanism. 
Furthermore, there are experimental indications that the 
symmetry of the order parameter of the high-Tc super- 
conductors could be an admixture of s and d- wave like 
(i.e, s ± d or d + is) Therefore, we present a model 
study for the coexistence of SDW with SC where the SC 
is caused by some generalised attractive potential with 
separable form (without specifying any particular mech- 
anism) leading to different order parameter symmetry. It 
therefore incorporates, two competing physical processes 
involving the electron-hole (SDW) like pairing of oppo- 
site spins with a net momentum difference {Q) between 
the conjugates and electron-electron (SC) pairing of op- 
posite spins with zero total momentum. Such a model 
study can also provide an important understanding re- 
garding the effect of the presence of antiferromagnetic 



(AF) order on the pairing symmetry of the supercon- 
ducting state, which is extremely important for high- Tc 
cuprates as well as the recently discovered borocarbides. 
Our main concern therefore boils down to study of the 
phase diagram, comprising of the gap magnitudes of the 
respective gaps e.g, the SDW and the SC gaps as a func- 
tion of band filling (hole concentration). Interestingly 
enough, we obtain the phase diagram which resembles to 
that of the high- Tc superconductors only when the super- 
conducting order parameter has the d-wave symmetry. 

The rest of the paper is arranged as follows. In sec- 
tion II we provide a brief mean field theory of the SDW 
ground state from a single band Hubbard model. Sec- 
tion HI contains the mean field model Hamiltonian for 
the coexistence of the SDW and SC together with the 
derivation of the respective gap equations. Section IV is 
devoted to discuss the results obtained from our model 
calculatons and finally we conclude in section V. 

II. SPIN DENSITY WAVE STATE 

The SDW state is a kind of antiferromagnetic state 
with the electronic spin density forming a static wave. 
The density varies perpendicularly as a function of po- 
sition with no net magnetisation in the entire volume. 
The SDW transition occurs when the spatial spin density 
modulation is due to delocalization or initinerant elec- 
trons rather than the localized one. Usually in the normal 
state the density p| (r) of electron spins polarized upward 
with respect to any quantization axis is completely can- 
celled by pi of downward polarized spins. In the SDW 
state, however, the difference (j{r) = P'\{'r)— Pi{r) is finite 
and modulate in space as a function of the position vec- 
tor r in the SDW state. Such tendency of forming SDW 
ground state takes place when a system possesses nested 
pieces of FS together with intermidiatc coulomb correla- 
tion. In the case of the SDW transition it is the wave vec- 
tor dependent static magnetic susceptibility which devel- 
ops a singularity at g = Q i.e., x(9,w)|^= Q j^^o = << 
S+{q,t) ; S~{-q,t) »^ |g=Q,a)=o — > oo where S"^ are 
the spin raising and lowering operators and Q is known as 
the nesting wave vector which determines the periodicity 
of the SDW. This singularity in the magnetic susceptibil- 
ity is an artifact of the nesting property of the FS given 
by 

= - 'k+Q (1) 

In addressing the possibilty of the coexistence of SDW 
and SC ordering in Cu02 based high- Tc material we will 
look for a model of a single Cu02 layer and neglect the 
inter-layer hopping because inter-layer hopping is weak 
compared to the intra-layer hopping integral. The simpli- 
est model one can use to describe the antiferromagnetic 
ordering of the two dimensional correlated system is the 
Hubbard model on a square lattice. 
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(2) where Asdw = —U ^ < C\.^Ck-Q[ >, the order pa- 



<ij>(y 



where tij is hopping integral between the different or- 
bitals ihj, fiia = Cj^Cja = number operator and U is 
the intratomic coulomb repulsion or the cost in cnc;rgy to 
put two electrons on the same site. The interaction term 
in the Hubbard Hamiltonian may be treated in the mean 
field approximation to get the SDW state as follows, 

Hi = UJ2ni^riii = uY^Cj^Ci^Cldi 

i i 
i 

- < C^Ci^ > cj^Cii- < cj^Cii > ClQ^] (3a) 



(3b) 



where Ui = i (< Cj^Cif > + < C^^C^ >), 
Af = -f(<4CiT> - <Cla^>)^^U<S.>, 



A+ = - U < C^C^ > = - U < S- > and 
A- = - U < Cla^ > = - U <S-> 

Therefore, the meanfield Hubbard Hamiltonian can be 
written as, 

H = -Y^ tij ClCj^ + UY1 ^i^i + (4) 



where hi = Cl^C,^+Cl^d and Bl = A| = - U < 
Sf > = the order parameter for longitudinal magne- 

transverse 



tization, 



1 , where 



± - A± = -V <Sf> 

magnetization and ct^ = ^ C\y ^ f 

f = Pauli matrices. The 2"'' term of equation (4) 
corresponds to total charge of the system whereas the 
3'''^ term is the total spin. Hence, eqn. (4) demonstrates 
that the charge and spin degree of freedom of the Hub- 
bard model is seperated. Now, the SDW state can be 
described equivalently in terms of either the longitudinal 
or transverse spin polarization. For example the SDW 
state with wave vector Q, having longitudinal spin po- 
larization is 

< 9^ >= X / cos((3..Ri) for SD on sites. , , 
^ "1 sm.{Q.Ri) for SD on bonds. ^ ' 

In momentum representation the mean field Hamilto- 
nian for a transverse SDW state becomes 

HSDW = '^i^k - f^)Cl,a^k,a + AsjjW X](*^fc-M3T^^i 



-h.c) 



(6) 



rameter for the transverse SDW state. The quasi parti- 
cle energy spectrum of the SDW state can be writen as 
Ek = \J {ek - m)^ + ^SDW "^ith a gap of 2 Asdw at the 
Fermi level. To note, there is a formal similarity of the 
SDW mean field theory with that of the BCS, in energy 
spectrum, gap equation, as well as in the collective modes 
(it is 2 Asdw in case of the SDW state). 



III. THE COEXISTENCE OF SDW AND SC 

Usually, magnetism and superconductivity are ex- 
pected to be mutually exclusive phenomena i.e, they 
are unlikely to occur simultaneously in the same com- 
pound. Superconductivity (including in the high temper- 
ature superconductors) is known to be due to Cooper pair 
formation of electrons of opposite spins and momenta 
whereas magnetism rcqiiircs spin polarisation. Therefore, 
naturally one order would inhibit the other. Further- 
more, like superconductivity (electron-electron pairing), 
the transverse SDW state is also a result of condensation 
of electron-hole pairs of opposite spins but with a mo- 
mentum difference of Q between the conjugates. Hence, 
when any of the orders (either the SC or the SDW) set 
in, the FS is instable with respect to that 'condensate 
state'. In other words, if one of the phases (say SDW) 
sets in first, and exists all over the FS, then there will be 
no carrier available to form cooper pairs and hence no su- 
perconductivity. This would also be equally true in case 
of the SDW state, had the superconductivity appeared 
first. However, in reality we do see the coexistence of 
the two phases as is already discussed earlier. In what 
follows we discuss two different scenarios where both the 
presence of SDW state and SC coexist. 



A. The case of an isotropic SDW-state : SC can 
arise over the SDW backg round 



If the nesting of the FS is perfect the entire FS could 
be isotropically gapped due to the formation of the SDW 
state transforming the system from a metal to an in- 
sulator. In order to build superconductivity over such 
an insulating phase one needs to dope the system with 
charge carriers. In the event of doping the system with 
holes, or equivalently removal of electrons from the lower 
filled (valence) SDW band, there will be deviation from 
perfect nesting of the FS resulting in a local suppression 
of the SDW gap. This is so because, the gain in elec- 
tronic energy resulting from the formation of the SDW 
state is lowered due to removal of electrons from the va- 
lence band. This local suppression of the SDW gap acts 
as a potential well for the injected hole, in which it gets 
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self-trapped forming the so called spin bag. On creat- 
ing two holes it is energetically favourable for them to 
dig a deeper well and stay together provided the two 
holes have opposite spins to avoid Pauli exclusion prin- 
ciple. This is however nothing but a local Cooper pair 
and if such bags with two holes of opposite spins move 
coherently, the system will be superconducting. This is 
the essence of Schrieffer's spin bag model for high tem- 
perature superconductivity jl2j and hence is an example 
where superconductivitycan arise over the SDW state. 



B. Anisotropic SDW state : SC can coexist with 
SDW 

The above case is however, physically resonable only 
when the hole concentration is very small. But in reality 
superconductivity in most of the systems (specially high 
Tc systems) appears only after large doping. So, it is 
likely that the SDW state will be completely suppressed 
in particular directions of the FS whereas it would still 
exist in the rest of the FS which still nests (anisotropic 
SDW). Such situation may also appear due to particular 
topology of the FS of a system without doping i.e, the 
system may have nested pieces of FS only in certain di- 
rection and no nesting in other direction resulting in an 
anisotropic SDW state. In the regions of the FS where 
the SDW gap vanishes (i.e, the lower and upper SDW 
bands merge together), the pairing interaction between 
the SDW quasi particles can take place leading to su- 
perconductivity. From a microscopic point of view, the 
effective attractive interaction between the SDW quasi 
particle in this picture can be rationalised (Ghosh and 
Sardar [|l^) as arising due to the exchange of the collec- 
tive modes of the SDW state. In both the cases discussed 
above, the origin of superconductivity is fundamentally 
different from that in conventional superconductors. 

In contrast, in the cases where nesting is not perfect 
either due to peculiar topology of the FS or due to dop- 
ing, the SDW gap will appear only in the nested part 
of the FS allowing for the superconducting instability 
in the rest of the FS. But in such a case, the origin of 
superconducting pairing may arise be due to any other 
mechanism including the BCS phonon exchange [p^ . Be- 
sides, at least in the case of high- Tc SCs there are other 
proposals for SC like the charge transfer mechanism fl^ ] 
and the pair tunneling mechanism etc. However, the 
nature of SC pairing and the symmetry of the order pa- 
rameter in the cuprates is still unknown. It was proved 
by Dzyaloshinskii and Yakovenko from a generalised 
interaction through renormalisation group equations that 
the two dimensional Hubbard model could lead to vari- 
ous instabilities including single superconducting dx2_y2 
pairing, along with possible coherent mixture of charge 
density wave (CDW) and SDW phases. Ruvalds and co- 
works 1 19 1 also found that (1^2 pa iring was favoured 



for nested fermi surfaces consisting of parallel orbit seg- 
ments. 

The Hamiltonian for the coexistence phase of super- 
conductivity and the SDW state can therefore be ob- 
tained by adding a pairing interaction term to equation 
(6) which in the mean field BCS approximation can be 
written as 



k.a k 

-Y,^sc{k){Clcl,^+h.c.) (7) 



where A^cik) 



k,k' 



Vk,k' < Cl^C^k'i >, Vk,k' being 



the strength of the attractive pairing interaction medi- 
ated by some boson exchange. Depending on the nature 
of pairing interaction the superconducting gap function 
exhibits different types of symmetry structure. 

In a weak coupling pairing theory, the SC gap param- 
eter satisfies the equation. 



= (^')tanh^ 



2Ek' 



(8) 



where = -^/e^ + A'^^{k) is the energy of the supercon- 
ducting quasiparticles. While Asc(fc) is not gauge invari- 
ant and therefore cannot be directly observed, j Asc(fc) p 
is observable. For a nodeless Afc(= Asc{k)), the state 
is non-degenarate and | A^ p has the symmetry of the 
crystal. However, if A^, has nodes in certain directions 
of k, there may be several degenerate pairing states and 
Afc p need not exhibit the point symmetry of the crystal 
— which is usually referred to as unconvensional pairing 
state. More precisely, the symmetry of the order param- 
eter is determined by the form of the potential V and the 
band structure. In the simpliest assumption the band en- 
ergies are independent of the direction in the a — 6 plane 
and hence the interaction potential can be approximated 
as a separable one, i.e, 



Vk,k' 



-Vf]kr]k' 



(9) 



where (a) rjk = constant, corresponds to an isotropic 
s-wave (BCS) symmetry ; (b) rjk = f(k), refers to an 
anisotropic s-wave symmetry, provided f(k) is a smooth 
function in the first Brillouin zone (BZ) and is -l-ve def- 
inite (i.e, nodeless) (c) rjk — cos(fc^a) -I- cos{kya), cor- 
reponds to an extended s-wave (s*) pairing symmetry 
whereas (d) rjk = cos{kxa) — cos(fcj,a) corresponds to 
dr^2_y2 pairing symmetry which also changes sign (like 
(c)) but does not transform as the identity under the full 
crystal group. 

Now, in order to study the interplay between the SDW 
and the superconducting state (with different possible 
pairing symmetry as discussed above) we need to calcu- 
late the self-consistent gap equations for the respective 
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gaps (e.g, SDW & SC). For this purpose it is convenient 
to write the Hamiltonian (7) in diagonahsed form using 
four component Nambu operator [|T] as follows, 

Hdiag = ^i'li'^kPscr^ + AscPaCTiAsDwPio-aHk (10) 



where the four component Nambu operators are defined 
as, 



and cr, p are 4x4 matrices defined as, 



(11) 



n 

T, 



(12) 



Ti being the 2x2 Pauli matrices and pi are given by 



Pi 



P3 
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(13) 



The gap equations corresponding to the order param- 
eters AsDW and A^c defined in equations (6,7) can be 
calculated by using equation (10) which are given below 
along with the equation for the total number of charge 
carriers. 



A 



SDW 



(4) E (-l)^^tanh.^^^(^) 



k,i=1.2 



and 



^^(^')tanh^ 



fe'i=1.2 



71 = 1 - 



- E 

2 ^ 

fe,i=1.2 



(e, 



'^^tanh^ 



(14) 



(15) 



(16) 



with 



E,ik) = y^(efe-Ai)2 + Af(fc) and 

A,;(fc) = (A,,(fc) - i-iyAsDw)- (17) 

The equations (14-16) are coupled integral equations 
in the sense that A^^ = ^sd^sc, ^soWjIt-), Asdw = 
^SDwi^sc^SDWin) and so is n = n{Asc, AsDW,n). 
The self-consistent solutions of these three equations are 
obtained numerically. Before the numerical results are 
presented it is worth pointing out that it appears from 
equations (14-17) as if there exists two effective order 



parameters Ai^2 = Age ± Asdw indicating that these 
can interfere with each other either destructively or con- 
structively. In striking contrast, in the coexistent phase 
of the charge density wave (CDW) and superconduc- 
tivity no show such interference is obtained where 
the quasi particle energy spectrum has an effective gap 
^ ^QY)^ + A^ , AcDW being the CDW order param- 
eter. Such a difference in behaviour between the SDW 
and CDW coexistent with SC may be attributed to the 
effect of interference between the order parameters of the 
SDW and SC states both of which involve the pairing of 
up and down spin quasi-particles, i.e, the electron-hole 
and electron-electron pairings respectively. Therefore, 
the electrons with the same spins are likely to compete 
for both the processes, thereby giving rise to interference. 



IV. RESULTS AND DISCUSSION 

The interplay between the SC and SDW state is mainly 
studied by selfconsistently solving the gap equations (14 
- 16) with fixed set of parameters C/ = 1, V ^ 1.5, the 
cut-off energy for SC = 0.8, all are in units of the 
hopping intgral t. Two distinct set of results are ob- 
tained, (i) for a fixed temperature (5 K) the amplitudes 
of the gap functions (SDW & SC) are obtained for dif- 
ferent pairing symmetry as a function of band filling (n) 
and (ii) temperature variation of the respective order pa- 
rameters (for different pairing symmetry) are obtained 
for different (fixed) band fillings. In the former case, the 
thermal variation of the chemaical potential, /i(T) is not 
important whereas for the later the thermal variation of 
p(T) is taken care such that given value of n differs from 
that of the calculated ones at best at the fifth decimal 
point. Furthermore, in computing the self-consistent so- 
lutions of equations (14 - 16) for different pairing symme- 
tries of the superconducting state, specific forms of the 
pairing potential and the corresponding order parameters 
are also assumed. Following the form of the potential in 
equation (15) we assume the form of the order parame- 
ters, for example, for the extended s (s*) and d wave SC 
as, A-'^'') = A°(cosA;j; ±cosky). 

Now we will discuss the results of our calculation. In 
fig. 1. we have plotted the isotopic s-wave SDW gap 
parameter (^Asdw^ and isotopic s-wave SC gap param- 
eter [Asc] evaluated at T = 5K vs. hole concentration 
X = \ — n. Clearly the figure demonstrates the strong 
competition between the two orderings namely SDW and 
SC emphasizing the coexistence of SDW state and Super- 
conductivity at low hole concentration x < 0.23 . As we 
increase the hole concentration first SC gap and then the 
SDW gap goes to zero. As a result 0.23 < a; < 0.25 wc ob- 
serve a pure SDW phase. When the hole concentration 
is further increased we obtained pure superconducting 
phase for 0.25 < a; < 0.5. This reentrance of supercon- 
ducting phase with respect to hole concentration is quite 
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a new phenomena and to the best of our knowledge no 
such theoretical demonstartion exists. 

In fig. 2. we have shown the variation of SDW gap 
function and amplitude of extended s-wave (s*) SC gap 
function A^^ evaluated at 5K with the hole concentra- 
tion X for the same set of parameters as in Fig. 1. Here 
there is no such coexistence of SDW and SC phase. In 
fact SDW phase is dominant in the low hole concentra- 
tion region where as SC phase is observed at large hole 
concentration. More precisely SDW gap shows a max- 
ima at half filling {x = 0), decreases with the increase 
of x and ultimately vanishes at around x=0.5 where as 
extended s-wave SC gap is seen at around x > 0.65 with 
a maxima around x = 0.8. Comparing the figure 1 and 
2 we find that superconducting ordering suppresses the 
SDW ordering and vice versa although for isotopic s-wave 
pairing one finds the coexistence of SDW and SC phase. 

In fig. 3. we have presented the scenario of d-wave 
superconductivity with the SDW. In the undoped (x=0) 
state the system is an insulating antiferromagnet with a 
SDW gap. As the system is doped with holes the an- 
tiferromagnetic ordering is suppressed and the system 
becomes superconducting in the region 0.4 > a; > 0.25. 
When X > 0.4 superconductiv ity disappears and we ob- 
tain a reentrant SDW phase. However, in this region of 
hole concentration {x > 0.4) the SDW gap parameter is 
small. Possibly the quantum fluctuation which has been 
neglected in this mean field calculation will destroy the 
SDW ordering. 

The SDW gap being isotropic s wave it has always 
maximum value at the half-filling because the density of 
states (DOS) is maximum at the Fermi lavel (half-filing) 
only. Since the Fermi surface is determined by the zero 
energy contour of ek = —2t{coskxa + cos kyu) = 0, the 
factor (cosfca;a-|-cosA;j,a) is therefore smaller, for a square 
lattice, close to half-filling (as n — 1 corresponds to the 
position of the Fermi Surface) and has larger value away 
from half-filling. Therefore, the self consistent solutions 
of the amplitudes of the extended s (d) wave supercon- 
ducting gap will always be minimum (maximum) at half- 
filling whereas maximum (minimum) at zero filling. This 
would qualitatively explain the band filling dependence 
of the extended s wave and d wave superconducting gaps 
(as is presented in Figures 1-3). The presence of the 
SDW order makes the situation complicated (and proba- 
bly more realistic for the systems like the cuprates) . The 
magnitudes of the SDW gap is maximum at half-filling 
not only because of the large van Hove DOS but also due 
to the fact that the nesting of the FS is strong close to 
half-filling. This results in the appearence of the SDW 
phase at lower doping concentrations (x) irrespective of 
the nature of the superconducting pairing (cf figures 1-4) . 
Considering the case of interplay between the SDW and 
s wave superconductivity we therefore find that both the 
SDW & SC states are degenerate at x = 0. However, 
the SDW phase is more stable than the SC phase al- 



though the interaction strength is smaller {U < V) due 
to strong nesting effect. Contribution to superconduct- 
ing phase at half-filling is probably from the non-nested 
parts of the FS. However, with doping the nesting of the 
FS gets affected and superconducting phase becomes sta- 
ble suppressing the SDW phase (cf Fig. 1, Fig. 3). In 
case of extended s wave pairing therefore, no strong in- 
terplay is found as expected (cf. Fig. 2). In case of d- 
wave pairing (cf. Fig. 3) however, the amplitude of the 
SDW order parameter is large at x = compared to that 
for s wave pairing scenario. A closer comparision of the 
figures 1 & 3 will also reveal that the maximum value of 
d wave SC gap is approximately three times larger than 
the maximum value of isotropic s wave SC gap (within 
the same set of all other parameters). This is because for 
isotropic s SC wave pairing, both the pairing strengths 
(SC as well as the SDW) are onsite in nature and will 
be most effective only at the half-filling due to symme- 
try reason (as discussed earlier). As a result both the 
orders will inhibit each other. In contrast to isotropic s 
wave scenario, for d wave pairing the SC phase is how- 
ever pushed towards higher doped region (unlike being 
degenerate at x = 0). This is due to the fact that the 
more extended d-wave pair wave function is more effi- 
cient in avoiding the pair breaking effect of the staggered 
local field. In the conserving fluctuation exchange ap- 
proximation (FLEX) which are in good agreement 
with available monte carlo data at higher temperatures 
showed that when the two dimensional Hubbard model 
was doped beyond a critical doping the leading instability 
changes from the SDW channel to d^^-y^ pairing channel. 
Therefore, our results in the flgure 3 from a mean fleld 
calculation agrees qualitatively with that of FLEX pC| ] 
approximation calculation. Furthermore we observe that 
the gap function for d-wave superconductors decreases 
sufficiently rapidly from its maximum value with hole 
concentration in comparison with isotopic s-wave or ex- 
tended s-wave superconductivity. Finally, in contrast to 
s-wave pairing d-wave superconductivity does not show 
any coexistence of SDW and SC phase. 

In figure 4. we have plotted the SDW gap and an ad- 
mixture of s and d wave pairing (s-fd) superconductivity 
(with SC gap having 80% d-wave contribution and 20% 
extended s-wave contribution) vs x. Here we find that 
the superconducting phase appears only within very nar- 
row range of hole (almost does not appear) concentration. 
Also the highest value of the SC gap amplitude is much 
smaller than that seen for the pure s or d wave picture. 
This is probably an indication of the fact that the su- 
perconducting state with mixed s + d symmetry of the 
order parameter is stable only in presence of orthrhombic 
distortion not in tetragonal phase ||2^. We shall present 
results in details of the interplay of the SDW state with 
complex (s + id) SC pairing in our next publications. 

Fig. 5 & 6 demonstrates the thermal variation of the 
s-wave SC gap and the SDW gap for diflerent hole con- 
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centrations. In Fig. 5, (for given set of parameters as 
earlier) for the half-filled case {x = 0.0) the SC gap opens 
up at higher temperatures {Tc =161K), grows with the 
lowering in temperature but drops to zero at lOOK till 
75 K. On the other hand, the SDW gap opens up at 
lOOK and grows with the lowering in temperature but its 
growth being arrested at 75 K (due to the re-entrance of 
the SC state at the same temperature) on further lower- 
ing the temperature. Away from half-filling, at a; = 0.2, 
the SC Tc is reduced drastically (to ~ 85K) and the SC- 
gap drops to zero below Tc at 28 K and reappears at 
21 K. The SDW gap (at x = 0.2) grows hke a first or- 
der transition at around 28 K (where SC-gap drops to 
zero) and grows with lowering in temperature which is 
again being arrested with the reentrance of SC below 21 
K. Therefore, Fig. 5 demonstrates very strong interplay 
and co-existence between the SDW and the s-wave state. 
To be noted the strong interplay between the two or- 
dered phases that includes their coexistence as well as 
re-entrance, is maximum at half-filling. 

While in Fig. 5 there exists coexistence as well as 
re-entrance of different parameters. Fig. 6 demonstrates 
another particular aspects of their interplay. The hole 
concentration here is x = 0.1 and one finds that for cer- 
tain range of thermal region, around 48 to 65 K none 
of the order parameters are stable. In other words, 
within this thermal regime both the order parameters 
precipetously either goes to zero or a finite value alterna- 
tively for a very small change in temperature. In the rest 
of the thermal regime it retains the essential features of 
the earlier figure (Fig. 5) as mentioned earlier. The arrest 
of the growth of the SDW gap at lower temperature with 
the re-entrance of the SC gap can be understood as fol- 
lows. The SDW gap (in the mean field treatment) being 
a very sensitive function of the density of state available 
at the FS, the opening up of the SC gap takes away large 
density of states and hence the SDW gap is suppressed 
at lower temperatures. On the other hand, the rentrance 
as well as the thermally fluctuating natures are due to 
typical interefernce between the two order parameters 
which are described to some extent analytically by the 
equation (17). Therefore, the present paper presents as 
extensive study on the influence of the SDW order on the 
SC pairing symmetry. We would like to further point out 
that the figures 1 & 5 have resemblence with the spin 
susceptibility data of borocarbides indicating an 

interplay of magnetism and superconductivity in these 
systems. 
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FIG. 1. Amplitudes of the SDW and the s-wave SC gap 
evaluated at 5K, as a function of the hole concentration (x). 



FIG. 2. The SDW and extended ,s-wave gap amplitudes 
evaluated at 5 K is plotted as a function of the hole concen- 
tration (x). 



FIG. 3. Gap amplitudes of the SDW and the d^2_y2 SG 
gap a s a function of hole concentration (a;). 



FIG. 4. The SDW and SC gap in the mixed (s + d) sym- 
metry is plotted as a function of the hole concentration. 



FIG. 5. Thermal variation of the isotropic s-wave SC gap 
and the SDW gap (in units of t) for various fixed hole con- 
centrations (x). 



FIG. 6. Temperature variation of the isotropic s-wave SC 
gap and the SDW gap for x = 0.1. The thermally instable 
region is worth noting. 
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